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Abstract 

The equivalence of the covariant renormalization and the partial- wave renor- 
malization (PWR) approach is proven explicitly for the one-loop self-energy 
correction (SE) of a bound electron state in the presence of external pertur- 
bation potentials. No spurious correction terms to the noncovariant PWR 
scheme are generated for Coulomb-type screening potentials and for external 
magnetic fields. It is shown that in numerical calculations of the SE with 
Coulombic perturbation potential spurious terms result from an improper 
treatment of the unphysical high-energy contribution. A method for perform- 
ing the PWR utilizing the relativistic B-spline approach for the construction 
of the Dirac spectrum in external magnetic fields is proposed. This method 
is applied for calculating QED corrections to the bound-electron ^-factor in 
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H-like ions. Within the level of accuracy of about 0.1% no spurious terms are 
generated in numerical calculations of the SE in magnetic fields. 
PACS numbers: 12.20.Ds 



I. INTRODUCTION 



The partial-wave renormalization (PWR) approach has been proposed a few years ago 
as a convenient but non-covariant method to perform the renormalization numerically 
in bound-state QED calculations. It has been successfully applied first in exact numerical 
calculations of the self-energy and vacuum-polarization correction of order a [[j],|],f|] and 
further applied in exact calculations of QED-corrections of order a 2 (a is the fine- 
structure constant). A fair agreement between the results obtained within different numerical 
approaches can be stated (see, e.g., the results for the effective self-energy correction in Refs. 
|||7]]). Nevertheless, questions about the equivalence between the covariant renormalization 
and the numerical PWR scheme respectively conjectures about the possible occurrence of 
spurious terms in numerical calculations of higher-order QED effects have been raised in 
the past H,f|. In Ref. || this issue has been anticipated qualitatively in connection with 
problems encountered in the numerical evaluation of the screened Lamb shift when non- 
covariant, numerical renormalization schemes are employed. Persson et al. || have made 
the first attempt to derive such spurious correction terms to the PWR analytically. To 
our knowledge this has been the first and only reference in which corresponding terms 
have been presented explicitly. They considered the exact self-energy correction of a bound 
electron state in the presence of an additional Coulomb-type screening potential V c which is 
treated perturbatively. Formulating the PWR by employing the Pauli-Villars regularization 
a generic, regulator-independent correction term that could contribute to the level shift of 
a bound state \a) is derived from corresponding counterterms (see Eq. (44) of Ref. ||): 




a 



(1) 
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Although for the particular situation under consideration the correction term cancels because 
it occurs with opposite sign in different subgroups of diagrams, the authors conjectured that 
this may not always be the case in calculations of higher-order QED effects and that the 
PWR and the covariant renormalization could lead to results that differ generically by 
terms of the form of Eq. (|IJ) . The conclusions drawn in || received further the support from 
numerical results for the SE including a perturbing 1/r potential ||10|| . 

In this paper we wish to address first the question about the occurrence of spurious 
terms of the generic type ([j]). Therefore, we reinvestigate the problem considered in Ref. 
||. In conclusion we find no indication for spurious terms generated by the numerical PWR 
method in the case of one- loop SE calculations in external fields. The correction (1) is shown 
to be due to an improper treatment of the unphysical high-energy contribution to the SE. 

Recalculating the examples, given in Ref. ]10| we find that the spurious terms originate 
again from a similar unphysical high-energy contribution. 

Finally, we investigate the problem of the spurious terms in an external magnetic field. 
In Ref. H it was conjectured that in the case of external magnetic perturbation spurious 
contribution to SE should remain. Employing a new approach developed here for the PWR 
in magnetic field, which is based on basis set expansion for the Dirac equation due to the 



Chen and Goldman |TTJ we prove the absence of the spurious terms for the QED corrections 
to the bound electron ^-factor numerically on the level of accuracy of about 0.1%. 



II. EQUIVALENCE BETWEEN THE COVARIANT AND THE PARTIAL- WAVE 

RENORMALIZATION 

In Ref. || a generic correction term between the PWR and the covariant renormalization 
contributing to the energy shift of a bound-electron state \a) interacting with a spherically 
symmetric perturbation potential V c has been derived within the Pauli-Villars regularization 
scheme. The authors obtain a spurious correction term to the PWR from both the A- 
dependent wave-function correction and from the vertex correction (see Eqs. (36) and (39) 
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of Ref. §): 

^ a (A) = -E/ ^ ^TTXT^ (aUl(fcr)(F e (r) -F c )|a) = — E / dk f^(k,A) ,(2) 

7T ^ JO k' (k + k'Y 7T ^ Jo 

where A denotes the Pauli-Villars regulator and fc' = V^ 2 + A 2 . Without going through 
the details of the derivations given in Ref. || we take Eq. (0) as the starting point of the 
following considerations. Note that the integral over the momentum k will be finite for each 
partial wave i. The k- and r-dependence of the integrand ensures a sufficient convergence of 
both integrals, which allows one to interchange the order of integrations. For any finite value 
r > the integral over k is sufficiently convergent, i.e., the integrand falls off as ~ 1/k 3 for 
asymptotic values of k. On the other hand the contribution to the matrix element arising 
from the integration over r from r = to some arbitrarily small value r = ro will be 
negligible. In the following considerations (part A) we keep the usual order of integrations 
as dictated by the PWR approach (see e.g. [0]), i.e. that the matrix element should be 
evaluated before the integration over k is performed. 

Suppose we could interchange the summation over i with all the integrations involved, 
then the correction term Eq. (0) vanishes evidently in view of the identity 

00 

E(2£+l)j 2 (A;r) = l (3) 
1=0 

after the integration over r is performed. Accordingly, the PWR approach and the covariant 
renormalization would be equivalent leading to identical results for the renormalized energy 
shift. Thus, one could try to prove explicitly whether or not Weierstrass' theorem for 
uniformly convergent functional series holds in case of the generic correction term (@). In 
the following we shall demonstrate that the functional series (|j) is uniformly convergent. 

A. Proof of uniform convergence 

The infinite summation over partial waves i may be decomposed into a finite sum < 
I < L — l, with L ^> 1 and the remaining infinite sum over I > L. Accordingly, it is sufficient 
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to focus on the remaining infinite sum. For this purpose we may substitute k = t/{\ — t), 
which transforms the indefinite integral over k involved in the generic expression (|2|) into a 
definite integral: 



In 00 r 1 
*l(A) = — E / * /?(*/(!-*), A) 

7T „_ r JO 



(4) 



As the next step we have to find an upper bound (majorante) u a t for each term ff of the 
functional series for all t e [0, 1] and for a fixed but large A>1 such that \ ff(t/(l —t),A)\ < 

wJ(A) and Y.T=l u e( A ) = C a ( A ) < 00 holds - 

If Weierstrass' criterium of uniform convergence is valid, we can interchange the summa- 
tion over i with the integration over t and thus over r as well. One key-point for finding a 
convergent majorante series is to estimate appropriately the square of the spherical Bessel 
functions involved. We rewrite Eq. (ffl) in the following form: 



££(A) —Y. I ' 



7T £ — JO 



,2+/3 {2£+l)(a\ 



\ rt ) JtKl-t, 



rP(V c {r)-V c )\a) 



(5) 



It 2 + (1 - t) 2 A 2 [t + ^t 2 + (1 - t) 2 A 2 
where (3 > 0. Note, that the integrand of the integral (5) has a complicated analytical 
structure when extended into the complex t-plane. Accordingly, the integral over t must 
be undestood as being performed along a suitable chosen contour from the very beginning. 
The choice of the contour of integration will be done below. The integral (5) is difficult to 
handle analytically in closed form. Therefore, we may now employ the following estimates 
for the two factors involved in the integrand: 



f 2 



t) 2 A 2 



t + Jt 2 



■i 2 



t) 2 A s 



< 



1 



A> 1 



(6) 



and for all £ > L ^> 1 and for all complex arguments z inside of a bounded region of the 
complex t-plane (see appendix A) 



\Ft(z,/3)\ 



+ 



[l-0(r e )], e>0 . (7) 



Here a' e l denotes the position of the first maximum of the spherical Bessel function jg and 
cp denotes some finite constant. We should point out that the latter inequality (0) strongly 



overestimates the function Ft by a constant value, which is even larger than the first (the 
largest) maximum of the function itself. Although the integral over t in Eq. (5) is convergent, 
one should be careful when treating the high-energy (t — *■ 1) region, which can actually 
now generate unphysical spurious terms as a consequences of the approximations performed 
above. We could restore or simulate the asymtotics by introducing some appropriately 
chosen convergence factor, e.g., e _Ati / < - 1 ~*' ) with /i > 0, if needed. 

Taking the absolute value of Eq. (||) together with the approximations (f|) and (|7|) we 
can write 



2a °° 
l^(A)l<-£ 

71 i=L 



1 t 2+/3 

dt — 
1 



t) 



1+/3 



Si 
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(a\r {V c {r)-V B )\a) 



(8) 



We observe the occurrence of Eulers' Beta function in the expression above as a consequence 
of the estimates performed. The Beta function is defined via the integral (see formulae (6.2.1) 
and (6.2.2) in Ref. 0) 



B(z,w) :-- 



dtt 



z-l 



t) 



w—1 



T(z)T(w) 



(9) 



T(z + w) 

Note, that for certain values of the arguments z and w the Integration over t has to be 
extended into the complex plane. Having in mind the analytical continuation, we can write: 



1 A ^ 



f 1 dt\B(3 + (3,-(3)\ = f 
Jo Jo 



i ^ r(3 + flr(i-fl 

PT(3) 



(10) 



The analytical continuation of the Beta function is provided by the reflection formula for 
the Gamma function (see (6.1.17) in Ref. [0): T(l — w) = —wT(—w) = ir/ sin(7rw) for 
< 3?{w} < 1. The rigorous treatment of the integral @, when evaluating it for the 
particular values z = 3 + (3 and w = —/3, is provided by contour integration (see e.g. Ref. 
|T3|). It is performed along Pochhammer's closed contour Cp on the Riemann surface of the 
integrand t z ~ l (l — t) w ~ l and relates the integral in Eq. (|8]), respectively, the Beta function 
(0) for arbitrary arguments z and w to a product of Gamma functions according to 



-wi (z+w) 



dtt 



z-l 



-ni (z+w) 



c P (1-t) 1 - 



J _|_ (z+w) g27i 



1 dtt 2 - 1 



— 4sin(7T2;) sin(7ru>) B(z, w) 



(11) 



Thus, we are lead to the final expression for the finite majorante series 

£ \ n A) I < |<V (K(r) - FJ| a >| J (/ + I)- 2 '" . (12) 

For 1/3 < (3 < 1 the sum over £ is convergent and can be expressed in terms of incomplete 
Zeta functions. Having derived a convergent majorante series the uniform convergence for 
the functional series (Q) is proven. This allows indeed the interchange of the summation 
over I with all the integrations involved. As result the generic correction term (0) is equal to 
zero. Thus, the spurious logarithmic terms derived by Persson et al. |§ do not occur. This 
demonstrates the equivalence between the covariant renormalization and the non-covariant 
PWR approach when applied to the problem under consideration. 



B. Comment on the missing term 

We would like to point out that the evaluation of the correction term (Q) as it has been 
performed in Ref. || is incomplete and that the logarithmic correction term Eq. ([!]) appears 
as an artifact of the way the expression has been evaluated. In Ref. |J, Eq. (0) has been 
evaluated according to limiting process: 



2a 00 

£?(A)+£ 2 a (A) = — £ lim 

7T e=Q ^oc 



Kr/r dk fl(k, A) + [ K dk fg(k, A) 

JKf/r 



(13) 



where f is some average value of the coordinate r. Following the arguments in Ref. || the 
second term vanishes in the limit K —>■ oo, while the first part £° generates the spurious 
term ([!]). In contrast to Ref. |J let us now take into account the third term, which is 
supposed to be zero when K tends to infinity: 

Consider the sum of the second and third term, i.e.: 



zm =— e i im / f^ A ) ( l4 ) 

7T „ n K^ooJk 



2rv 00 r°° 
£«(A)+£«(A) = — £ lim lim / dk e^ k ft(k, A) (15) 
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The regulator e~^ k is introduced for reasons of simplicity in order to guarantee a finite 
integral over k at the upper integration limit, if the factor jf is absent. For finite values 
of the parameter fi this regularization will generate some large but r- independent constant, 
which, however, cancels in the matrix element. Similarly, one may include the regulator 
e -/V fc j n t derive the same result (|I]). 

Now we can employ the same arguments that were used in Ref. || for calculating 
The authors of |§ argue that the order between the integration over r and the summation 
over I can be interchanged for the following reasons: a) the r-independent part of the k- 
integrand does not contribute to the matrix element (a\ ■ ■ ■ \ a), and b) the convergence of 
the k- integral then does not depend on the Bessel function jf . 

The rigorous treatment of interchanging summations and integrations requires the valid- 
ity of Weierstrass' criterion as we have demonstrated in part A. Employing the arguments 
a) and b) above, we assume that in Eq. (pT5|) the summation can be interchanged with all 
integrations and limits involved. Using the identity Eq. (|3|) we arrive at 

2a _ f°° A; e~ 5 



£2(A) + £«(A) = — lim lim(a|(K(r) - V e ) f 

7T if— >oo/i— >0 Jjt 



Kfii/r 



n 2 



| a) 



oi 

— lim lim (a \(V c (r) - V c 

Z7l K^oo fi^O 



r J \ r J 



V £ ) 
a) . (16) 



To obtain this result we performed the limit /i — > in the integrand first. Expansion of the 
exponential integral for small arguments: E\(^ ii ) = (—7 — In j + O r^jp) and taking 
the limit A — ► 00 of Eq. (|16|) one ends up with 

r 2 +£« = ^-(a\{V c {r)-V c )Hr)\a) . (17) 

Thus, we obtain the logarithmic potential term as the remaining cutoff-independent con- 
tribution. It carries an opposite overall sign and will cancel in the total sum. Moreover, 
we could introduce a unique regulator e _M< - 1//fc+fc ' ) in expression (Q) and thus avoiding any 
(unnecessary) decomposition of the /c-integral and the introduction of f. Interchanging and 
performing the summation over I at first one is left with 
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£ a (A) = lim (a\ (V c (r) - V c )2 \k,{2^) - A 2 K 2 (2fi) + ^A 4 K 4 (2/x) + • • •} |a) = . (18) 
Zir m->o ^ lb J 

All terms in the curly brackets vanish when the matrix element is evaluated. The calculation 
performed above indicates that the occurrence of the spurious logarithmic term (and may 
be others as well) strongly depends on the analytical and on the numerical treatment of the 
high-energy contribution to the SE. 



III. A NUMERICAL ANALYSIS ON SPURIOUS TERMS 



In Ref. |1(J the occurance of spurious terms has been reported in connection with nu- 
merical evaluations of the SE correction in one-electron ions with nuclear charge Z in the 
presence of additional Coulombic perturbation potentials. We will show, however, that this 
is due to an improper treatment of the high-energy contribution to the SE. The total ex- 
pression for the SE correction of a bound-electron state \a) can always be represented in the 



form (see |TJj] for details): 

AE^ E = ^(an|S|no), (19) 

n 

where £ denotes the nonlocal electron self-energy operator and the labels n run over the 
complete Dirac spectrum. The SE correction in the external field together with an additional 
perturbation potential can be divided into three parts : the wave-function correction, the 
vertex correction and the derivative (or reference state) correction j8|. We will concentrate 
here on the wave-function correction due to the self energy (SE, WF). This correction can 
be obtained from the lowest order AE^ E correction by a replacement of the unperturbed 
wave function \a) by its first-order perturbation theory correction 

i , v-^/ \n)(n\V\a) 

l«)->E E _ E . 

n - L -'a J - J n 

where V denotes the perturbation potential. The prime indicates that the term with E a = E n 
is omitted from the sum. The energy correction AE^ E,WF can be written as 
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AjBfl = ^ 1 E~^E + E < • (20) 



Within the B-spline approximation the complete Dirac spectrum is represented by a purely 
discrete one and terminates at some large number N = 4L(u + s — 2), which is determined 
by the number of partial waves L, the order s of the B-spline functions and the number of 
grid points v. We may restrict the summation over the energy of the intermediate states 
in both Eqs. (19) and (20) by the condition \E n \ < k ■ mc 2 . The numerical results for the 
energy corrections as a function of the cutoff parameter k are shown in Table 1. What the 
lowest-order SE (Eq. (19)) is concerned its exact value in U 91+ (Z=92) is already obtained 
for k = 5 within an accuracy of about 0.7%. Any further enlargement of the summation 
interval does not lead to any improvement of the accuracy. We should note that the accuracy 
also depends on the number of partial waves taken into account and on the number of grid 
points. In our case we typically used L=6, and v ~ 140 together with spline functions of 
order s = 9. The quoted accuracy is sufficient for our purposes since the contribution of the 



spurious term obtained in |T(J is supposed to be much larger (about 20%). 

From Table 1 we conclude that the summation interval with k = 5 for H-like U +91 
provides a sufficient approximation (about 2% deviation) for AE% ' . However, we have 
found that any further enlargement of this interval generates the spurious contribution as 
obtained in [TO]. From the physical point of view the 20% contribution that originates 
from the high energies larger then 5mc 2 (compared to the binding energy 0.3mc 2 ) is hardly 
understandable. Therefore, we conclude that the occurance of the spurious contribution in 
the numerical calculations reported in Ref. |TOJ has a similar origin as the spurious term (1) 
derived analytically in Ref. ||, i.e. an improper treatment of the high-energy contribution 
to the SE within the PWR approach. 
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IV. CALCULATION OF QED CORRECTIONS TO THE G-FACTOR. 



Now we show that the conclusion drawn in Ref. |§ concerning the inapplicability of the 
PWR approach to SE calculations in external magnetic fields also does not hold strictly. For 
this purpose we employ a new approach to the evaluation of the SE including an external 
magnetic field in the Dirac equation from the beginning. 

Accordingly, the problem reduces again to the evaluation of the lowest-order SE (Eq. 
(19)). The Dirac equation with an external magnetic field is solved by means of an approach 
due to Chen and Goldman which we discribe briefly in Appendix B. 

Within this approach the vacuum-polarization and self-energy corrections to the electron 
g-factor in hydrogenlike heavy ions will be calculated. At first we calculate the vacuum- 
polarization effect in order to determine the values of the magnetic field strength where this 
method remains stable. After this we turn to the calculation of the self-energy correction in 
the external magnetic field within the PWR method. 

A. Vacuum-polarization corrections to the electron ^-factor of hydrogenlike heavy 

ions 

Vacuum-polarization (VP) corrections to the electron ^-factor for H-like highly-charged 



(HCI) have been calculated within the the Uehling approximation For calculations 

of energy level of HCI this approximation is valid with an accuracy of about 10% for all 
Z values. Within the Uehling approximation the VP correction of a bound-electron state 

\a(B)) including the external magnetic field is determined by the matrix element: 

Z „, „ ,„„ / Z 



AE^(B) = (a(B)\ - -S(r)\a(B)) = (—S(r)) , (21) 

r V r J a(B)a(B) 



where 



SW^fe^f.^^^ (22) 



^ f°° e -2r X /a f x + M VX 2 + 1 

3n Ji \ x 2 J x 2 

and a denotes the fine structure constant (in atomic units a = 1/c) . Then the correction to 
the ^-factor results as 
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5g y P = L a ^_ a V ^ (23) 

where /xg is Bohr's magneton (in atomic units \ib = §)■ 

The results of the calculations for the ls-ground state of H-like ions in comparison with 
perturbation-theory results obtained within the Uehling approximation [|TB| are given in 
Table 2. An agreement better than 0.01% has been found for all Z values from 5 to 90. In 
these calculations values for the magnetic field strength B = \B\ within the range of 1.0 a.u. 
< B < 100.0 a.u. have been used. The numerical results are stable within this range of B 
values. 



B. Self-energy correction to the electron g- factor in hydrogenlike heavy ions 

For the evaluation of the self-energy correction in the external magnetic field we employ 
the original PWR scheme developed in |I],]2[] in combination with the approach due to Chen 
and Goldman |TT[ for solving the Dirac equation in cylindically-symmetric external magnetic 
fields. 

The formula for the SE correction in an external magnetic field for the state a(B) reads 



7T 



n(B) \ ri2 J ° K ~ En ( l ~ i0 ) + Ea J a(B)n(B)n(B)a(B) 



where the sum runs over the total Dirac spectrum, a denote the Dirac matrices, r 12 = |fi — r 2 \ 
and 5m a abbreviates the counterterm. According to [|T],0 this counterterm follows from Eq. 
(24) by replacing the summation over bound-states n(B) by a corresponding integration 
over free - electron Dirac states and by replacing the bound-state energy E n (B) in the 
denominator by the free-electron energy e. The bound states a(B) are expanded in free 
electron Dirac states \p,e > (where p is the electron momentum). The correction to the 
g-factor that arises from the SE contribution reads 

SE _ 2[Agf(g) - AEf(O)] 
d9a ~ Js~B • (25) 
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Numerical results for the groundstate of H-like ions with different Z values are given in 
Table 3 in comparison with data obtained via perturbation theory |16j and from the aZ- 
expansion (for low Z), respectively. As for the VP corrections numerically stable results are 
provided for values of magnetic field strengths within the same range 1.0 a.u. < B < 100.0 
a.u. For all Z values the deviations from the perturbation-theory results turn out to be 
smaller than 0.1%. 

Thus, contrary to the statement made in ||, this provides numerical evidence for the 
absence of the spurious terms in calculations of the SE correction in external magnetic fields 
performed within the PWR approach field at a level of accuracy better than 0.1%. 

We should note that the accuracy achieved in our approach will be not sufficient for 
obtaining accurate values for bound-state QED corrections to electron g-factors. The reason 
traces back to the fact that bound-state QED corrections are obtained via subtraction of the 
free-electron QED corrections from the values given in Tables 2 and 3. This leads to severe 
numerical cancellations that deminish the accuracy of the net result significantly. However, 
it does not effect our conclusions concerning the spurious terms. 



V. CONCLUSIONS 

Summarizing, we have provided evidence for the absence of the " renormalization cor- 
rections" for the particular case of the SE in external fields. Contrary to the statements 
made in ||, Jl0| we have found that via a proper treatment of the high-energy contribu- 



tion to the SE such corrections can be avoided numerically and analytically. However, the 
proof presented here cannot exclude the possible existence of spurious corrections in more 
complicated situations, e.g., in connection with high-order QED corrections. In the partic- 
ular case of the "loop after loop" second-order electron self-energy correction a discrepancy 
between results obtained within the PWR |I7j approach and from the covariant remormal- 
ization scheme |18| and [^] has been reported. This may indicate that this problem requires 



further investigation. 
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APPENDIX A: THE UPPER LIMIT FOR BESSEL FUNCTIONS. 



Let us consider at the beginning the real functions Fg(x,/3) = j|(x)/x /3 for values (3 > 
and i > L ^> 1 over the interval x G [0, a^i]. We adopt the notation of Ref. [ 12T| . Let a^i 
and a' t 1 denote the location of the first zero and of first maximum of the spherical Bessel 
function jg(x), respectively. We employ the formulae (10.1.59) and (10.1.61) of Ref. |12| : 



1 

+ 2 



In -2/3 

1 + 0.80SG(£+- 



1 - o[t 



jH^,i)~°- 8458 (^+2 



-5/6 



1 - 0(t 



(Al) 
(A2) 



with e > 0. Obviously, the zeros of the functions F% and je coincide. For any I > and 
(3 > the function Fe has exactly one maximum within the interval x G [0, a^i] located at 
x = b' ei , i.e., F' t {b' tl) (3) = and F"(b'g X ,(3) < 0. Derivatives with respect to x are indicated 
by primes. For the first derivative F[ evaluated at a' ei , we find 



'2xj e (x)j e (x) - (3jj(x) 



x 



13+1 



3tA a Li) 



-(3 K l ; x 1 < . (A3) 



This implies: 1. b' f 1 < a' e x and 2. Ft decreases monotonically over the interval b' tl < x < a' tl . 
Using the differential equation for the je we derive for the second derivative 

iV 



x? 



13+2 < U ' 



[a 



(A4) 



Llj 
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which holds for all values < (3 < L, since £ 3> 1. The fact that F"(a' e 1 , 0) < reveals that 
F( has its point of inflection somewhere between a' ei and a^i, i.e. Ft is a convex function 
over the interval Accordingly, we can estimate the value of Fi at the position b' ei 

of its first maximum, which finally provides an upper bound for the function Fe throughout 
the whole range x G [0, oo): 



F e (b' £>1 ,i3) < F e (a' itl ,l3) + 1*2(0^,0)1 • \b' iA - a' e J < (1 +/3) F e (a' iA ,f3) . 
Here we used Eq. ( |A3| ) and the fact that \b' ix — a' tl \ < a' tl . Finally we arrive at 

-5/3-/3 



(A5) 



F e (x,P)<(l + (3) J -f^-< C(i U + - 



i - o(r e ) 



(A6) 



What the proof is concerned we shall restrict to values < (3 < 1 only. 

To provide the function F as required along the Pochhammer's contour we need to extend 
the considerations above into the complex plane. The function F is analytic in the complex 
plane z = x + i(. A suitable region of the complex t-plane (see Eq. (4)) which includes 
Pochhammer's contour is depicted in Fig. 1. 

For small imaginary parts < e 1 all steps of the derivation go through for complex 
arguments z = x + i( with \(\ < e. The Taylor expansion yields 



\F{*,P)\ 



F(x,/3)+iC (J^F(z,{3) 



C=o 



+0(C 2 ) < \F{x,P)\+e F\x,(3) +0{e 2 ) , (A7) 



where < x < 1. Now we might define sup( 0<a . <oo \|F'(x,/3)| = M < oo together with 
e = (t + l/2) _5 ^ 3 ~^. In view of Eq. (A6) we can write 

(A8) 



, i x -5/3-/3 






l - o{r e ) 
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FIGURES 
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FIG. 1. The region of the complex t-plane for which the approximation Eq. (A6) is valid. It 
includes the Pochhammer's contour as depicted. 



APPENDIX B: FINITE BASIS SET SOLUTION OF THE DIRAC EQUATION 
FOR ATOMIC ELECTRONS IN EXTERNAL MAGNETIC FIELDS 

The Dirac equation for a bound electron in an additional external magnetic field reads 

= Em , (Bl) 

with 

H = H + H m (B2) 

and 

H = a-p + pm + V c (f), (B3) 

H m = ±a-[Bxr], (B4) 

where a, (3 denote the Dirac matrices, m is the electron mass, B is the magnetic field 

— * 

strength. The magnetic field is supposed to be directed along the z-axis: B = B e z . Vc(f) 
is the Coulomb potential of a (point-like or extended) nucleus . In Eqs. (B1)-(B4) atomic 
units are used. 
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According to Ref. [0] the variational solution of the Dirac equation (Bl) is obtained by 
means of trial functions 



^raax Kmax 

*"(f) = E E ^T(r) 

1=1 K 



(B5) 



The electron wave functions ^(r ) in the magnetic field possess cylindrical symmetry 
and can be expanded with respect to a finite basis set of the functions ^/^(r ) of spherical 
symmetry. The index k denotes the Dirac angular quantum number, /i corresponds to the 
total electron angular momentum projection, a K are the variational coefficients, and 2N 
defines the number of the basis-set functions. 

The next step employs the B-spline representation of the functions ^/^(r ): 



2JV 

E bh®J 
j=i 



(B6) 



with 



V 



o 



(B7) 



/ 



N+n\ 



o 



\ 
) 



:bs) 



y -i-X-^(r/r) 

where n = 1, . . . ,N and denote the B-spline representations of the finite basis set of 
radial functions |^0[ and Xn^^lr) are the usual spherical spinors. 

The variational solution of the Dirac equation (Bl) with the trial functions (B5) reduces 
to the diagonalization of the Hamiltonian (35) within the finite basis set defined by Eq. 
(B5)-(B8). As the result one obtains the full set of solutions of the Dirac equation for the 
atomic electron in an external magnetic field. 

In particular, the matrix elements of the operator H m with the wave functions (B6) - 
(B8) are given by: 

B 



H„ 



N+n 



■ r KK , A i , 

2 nn kk ' 



(B9) 



17 



where 

/ ^"OO 

rV£(r)V#(r)dr, (BIO) 



nn 







4k2 - 1 , ' " ' (Bll) 

sgn(^) ^gj^ ,for «' = -« 

It is convenient to rearrange the summations in Eq. (38) - (39) in such a way that 
the summation over partial waves L will be performened as the last one. In the B-spline 
approach one starts to count the states from the lowest negative energy state (J = 1). Then 
the state with J = N + 1 corresponds to the positive energy ground state lsi/ 2 . 

Within the B-spline approach the H-like ion is included inside a spherical box of radius 
-Rbox ~ Y a - U -- The number of grid points was N g = 150 and the order of splines k = 8. 
This corresponds to 2N = 2(N g + k — 2) = 314 energy levels that represent approximately 
the Dirac spectrum. With this choice the inaccuracy of the spline approximation for the 
ls\/2 state compared to the Chen and Goldman variational solution [llj becomes less than 



10~ 8 . To test the accuracy of our approach we have calculated the Zeeman splitting of the 
Is 1/2 state in the hydrogen atom for different field strengths B. 

In Table 4 the results for the corresponding Zeeman splittings are compared which have 
been evaluated within the approach Chen and Goldman and by means of perturbation theory. 
The latter have been obtained employing the standard formula |^TJ . The comparison reveals 



that for field strengths B up to 2 • 10 2 T the deviation from the perturbation-theory results 
is about 10~ 8 while for a field strength of about 2 • 10 4 T the deviation increase up to 10 -3 . 
The latter is due to the strong distortion of the atomic structure by the magnetic field. 
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TABLES 

TABLE I. AE% E and A£^ E ' WF for U 91+ as a function of the energy cutoff k of the negative 
and the positive Dirac spectrum (see Eqs. (19-20)). All energy values in a.u.. 





covariant 


k = 5 


k = 10 


k = oo 


PWR 


AE™ 


13.0714° 


12.9781 


13.2190 


13.2114 




A£ ;SE,WF 


0.4626 b 


0.4671 


0.4914 


0.5695 


0.5598 c 



a Ref . H| 
b Ref. |6| 
c Ref. p| 



TABLE II. VP corrections to the electron ^-factor for the ground state of hydrogenlike ions. 



$9i s i this work 



6g? s p , Ref. @ 



1 

5 

10 
20 
50 
90 



-3.0 TO" 11 
-4.0 -10~ 9 
-6.4-10" 8 
-9.37-10" 7 
-3.32-10" 5 
-4.3991-10- 4 



-4.2-10" 9 
-6.37-10~ 8 
-9.41-10" 7 
-3.32-10" 5 
-4.3995-10" 4 
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TABLE III. SE corrections to the electron g-factor for the ground state of H-like ions. 



Z 5gff, this work 5gff, Q dgff aZ-exp 

1 (2.31±0.01)-10- 3 2.322840-10" 3 2.32284-10~ 3 

5 (2.32±0.01)-10~ 3 2.323389-KT 3 2.323388-10~ 3 

10 (2.321±0.009)-10- 3 2.325472-10" 3 2.3249028-KT 3 

20 (2.332±0.007)-10~ 3 2.33692-10~ 3 

50 (2.469±0.012)-10~ 3 2.47162-10" 3 

90 (2.997±0.011)-10~ 3 3.04516-10" 3 



TABLE IV. Results for the Zeeman splitting AE% for the ground state in neutral hydrogen 
obtained in the Chen and Goldman approach (CGA) Jl ij and in perturbation-theory (PT). All 
values are given in atomic units (a.u). For the magnetic field strength: 1 a.u. = 2.35 • 10 5 T. 
B AE Z (CGA) AE Z (PT) 

0.5000066566 0.5000066566 
10" 5 0.5000116565 0.5000116565 
10~ 3 0.500506398 0.5005066477 
10" 1 0.5497433 0.5500057691 

1 0.872133 0.9999977813 
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